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Right CFD Model for the Right Turbulence Problem 


In an apocryphal story, Heisenberg said, “When I 
meet God, I’m going to ask him two questions: Why 


relativity? And why turbulence? I really believe he 
will have an answer for the first.” [1 


With this in mind, it’s not surprising that nothing in 
the computational fluid dynamics (CFD) world is 


more contested than turbulence models. 


It’s paramount that you have the right model, built 
for the right job. As a result, it’s important to know if 
the simulation tool you are purchasing has the 
models you need for your applications. This eBook 
aims to provide a good place to start. 


Regardless, finding the right model will take trial and 
error, and even some physical testing. Once you 

have your model you will discover the power of CFD 
to reduce development times and prototypes. 
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What to Expect in This eBook 


This eBook will: 


Laser-induced, fluorescent false-color image of a 
submerged turbulent jet. (Image courtesy of C. 
Fukushima and J. Westerweel, Technical University 
of Delft.) 


* Define turbulence model families, their applications and limitations 
¢ List turbulence models offered from the list of companies below. 


Note that other big names in CFD were asked to participate and either declined or failed to reply. 


What to Expect in This eBook 


Applications 
Internal flows and electronic cooling 


Aerodynamics (transonic flows) 


Models 


Wall treatment models 


RANS Models: Spalart-Allmaras (SA) 


General purpose, external flows with complex 
geometry 
General purpose, internal flows, jets, curvatures, 
separated flows 


Cyclones, swirling and other complex flows 


Flows in between laminar and turbulent 


Thermal fatigue, vibration, buoyant flows (ship 
design) 
External aerodynamics, aeroacoustics, bulk or wall 
turbulence 


Fundamental research and turbulence model 
creation 


Separation zones, aerospace, channel pipes 
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Wall Treatment Models 


¢ Specifications 


RANS Models: k-e 


RANS Models: k-o 


RANS Models: Reynolds Stress Transport 


Transition modeling 


LES Models 


DES, DDES, IDDES models 


Direct Numerical Simulation (DNS) 


Scale Adaptive Simulation (SAS) 


¢ Dimensionless wall distances (Y+) or length-velocities 


(L-Vel) calculate turbulent viscosity from local fluid 
velocities and the distance from the closest wall 


¢ Solves the flow everywhere 
¢ Very robust 
¢ Lowest computational cost 


¢ Applications 


* Internal flows 


¢ Particularly electronic cooling 


¢ Limitations 
¢ Typically has low accuracy 
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Wall Treatment Models 


ESI Dassault 
Models Altair Autodesk COMSOL ANSYS pan Systémes 
AcuSolve CFD Multiphysics SIMULIA 
OpenFOAM 
ABAQUS/CFD 
High Reynolds 
Number (High Y+) x x x 
Low Reynolds 
Number (Low Y+) x x 
All Y+ Hybrid 
7 z X X X 
reatment 
Length-Velocities xX 
(L-Vel) 
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RANS Models 


¢ Reynolds-averaged Navier-Stokes 


(RANS) is a large family of turbulence 
models 


¢ These models supplement the Navier- 


Stokes equations by adding viscosity 
terms 


¢ The added term is calculated 


differently for each model 


¢ The common term is k, the kinetic 


energy per unit mass of turbulent 
fluctuations 
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Comparison of results for the turbulent fluid flow 
through a prism. One is an instantaneous look and 
the other based on Reynolds-averaged Navier- 
Stokes equations. (Image courtesy of Emilio 
Baglietto.) 


RANS Models: Spalart-Allmaras 


¢ Specifications 
¢ One-equation turbulence model 
* No wall functions 
¢ New variable added to represent SA viscosity 
¢ Not memory-intensive 
¢ Stable with good convergence 
¢ Stanford notes that this is a good compromise for accuracy vs. 


simplicity 
¢ Applications: 


¢ Aerodynamics flows 
¢ Transonic flows over airfoils 


¢ Limitations: 


¢ Solving shear flows, separated flow, decaying turbulence 
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RANS Models: Spalart-Allmaras 


ESI Group / Dassault 

: rou : 

Models Altair Autodesk CoMSer ANSYS_ CED Direct Systémes 

AcuSolve CFD Multiphysics OpenFOAM SIMULIA 
ore ABAQUS/CFD 
Standard 
formulation 

(low-Reynolds x x x x x 
number) 

Rough-wall SA XxX 

High Reynolds X X X 


number 
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RANS Models: k-e 


¢ Specifications 
* Based on k and the dissipation rate of turbulent kinetic energy, € 
¢ Uses wall functions so buffer region flow isn’t simulated 
¢ Very popular two-equation model 
* Reliable, good convergence and low memory requirements 
¢ Numerous variants improve certain limitations of the model 
¢ Stanford prefers realizable k-e, Kato and Launder variations 


¢ Applications 
¢ General purpose: no separation, compressible/incompressible 
¢ Particularly external flow interactions with complex geometry 
¢ Limitations 


* Not accurate for no-slip walls, adverse pressure gradients, strong 
curvature into flow, and jet flows. ¢ is difficult to solve. 
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RANS Models: k-e 


ESI Group / Dassault 
Models Altair Autodesk COMSOL ANSYS CFD kee Systémes 
AcuSolve CFD Multiphysics OpenFOAM SIMULIA 
Ber ABAQUS/CED 
Standard k-e xX xX xX x 
Standard two- 
x 
ayer k-¢€ 
Standard k-e low- 
xX x 


Reynolds number 


Scalable Wall 
Function k-¢ x x 


Realizable k-e xX xX xX x x 


Realizable two- xX 
layer k-e 
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RANS Models: k-e 


ESI Group / Dassault 
Model Altair Autodesk COMSOL ANSYS CFD ses Systémes 
va AcuSolve CFD —— Multiphysics OpenFoAM _SIMULIA 
Be ABAQUS/CFD 
RNG two-layer 
7 xX 
-€ 
RNGke X xX xX xX x 
Lien Cubick-e  X 
V2F low-Reynolds Xx 
number 
SSG Reynolds 
stressed k-¢ x x 
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RANS Models: k-¢ 
Dassault 
ESI Group / 
Altair Autodesk COMSOL : Systémes 
ANSY FD D t 
maouels AcuSolve CFD Multiphysics “NS¥S . ona SIMULIA 
vs ABAQUS/CED 
Rough-wall k-¢ xX 
xX x 


LRR k-e 


Launder and 
Sharma Rapid xX 
Distortion Theory 


GibRon ner’ 
Dafa’ Alla’s q-zeta 


two-equation xX 
low-Reynolds 
number 
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RANS Models: k-e 
D It 
Altair Autodesk = COMSOL ESIGroup/ tam 
air utodes : ystémes 
Model ANSY FD Direct 
“_ AcuSolve CED ‘Multiphysies. “™°** 2 a SIMULIA 
pen ABAQUS/CFD 
Bouyant k-e xX 
Multiphase x x 
standard k-e 
Multiphase 
standard two- xX X 
layer k-e 
Multiphase xX 
Lahey k-e 
Multiphase xX xX 


realizable k-¢ 
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RANS Models: k-@ 


¢ Specifications 
¢ Based on k and the specific rate of dissipation of kinetic 
energy, @ 
* © used as it is easier to solve than € 
¢ Uses wall functions, so buffer region flow isn’t simulated 
¢ Similar memory requirements as k-g 
¢ Stanford suggests that the STT variation of k- @ is 


typically better than the standard 
¢ Applications 


¢ Similar to k-e with improved accuracy for: 
¢ Internal flows, curvatures, separated flows and jets 


¢ Limitations 
¢ Hard to converge vs. k-g, sensitive to initial conditions 
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RANS Models: k-@ 


ESI Group / Dassault 
Altair Autodesk COMSOL : Systémes 
Y FD Direct 
eee Acusalve’ CRD ‘Muliphysies “°°” : 7 Aa SIMULIA 
ove ABAQUS/CFD 
Standard k-w xX xX x xX xX 
SSTko X xX X xX xX XxX 
SST Rough-wall p 4 
Baseline (BSL) 
k-@ x x 
Rough Wall k-o xX 
Multiphase Sato xX xX 


SST k-@ 
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RANS Models: Reynolds Stress Transport 


¢ Specifications 
¢ Attempts to directly model flow and terms in RANS equation 
¢ Based on six equations for turbulent stresses 
¢ Good physical representation of flow 


¢ Applications 


¢ Good for assessing new phenomena and complex flows 


(cyclones and swirling) 


¢ Limitations 
¢ Computationally expensive, sensitive to initial conditions, lots 


of modeling, needs high-quality mesh 
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RANS Models: Reynolds Stress Transport 
Dassault 
ESI Group / 
Altair Autodesk COMSOL . Systémes 
ANSYS_ CFD Direct SIMULIA 


mode AcuSolve CFD — Multiphysics OpenFOAM 
Pee ABAQUS/CFD 


XxX 


Linear 
pressure strain 


Quadratic 
pressure strain 


Linear pressure 
strain two-layer 


Explicit Algebraic 
Reynolds Stress xX xX 
Model (EARSM) 


Low Reynolds 
EB-RSM 
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Transition Modeling 


¢ Specifications 


¢ Models how the flow transitions from laminar to turbulent 
¢ Transition isn’t discontinuous 


¢ Dynamic process that is situationally dependant 


temporally, spatially and spatio-temporal 


¢ Applications 


¢ Flows that move between laminar to turbulent 


¢ Limitations 
¢ Modeling and documenting boundary conditions making 
acoustics and disturbance modeling difficult 


Page 20 


Transition Modeling 


Dassault 

ESI Group / 
Altai Aut k é 
air utodes COMSOL ANSYS_ CED Direct Systémes 


AcuSolve CFD Multiphysics SIMULIA 


OpenFOAM 3 AQUS/CEFD 


Models 


User specified 
turbulent xX 

suppression 

Gamma ReTheta 
X X 

transition model 
SST + 1 equation 
transition (SST- xX 
gamma) 
Low Reynolds 
Number k-kL k-@ x x 
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Transition Modeling 


ESI Group / Dassault 
Models Altair Autodesk COMSOL ANSYS CED Direct Systémes 
AcuSolve CFD Multiphysics OpenFOAM SIMULIA 
ABAQUS/CFD 


Spalart-Allmaras 
+ | equation X 
transition model 
(SA-gamma) 
Spalart-Allmaras 
+ 2 equation 
transition model xX 
(SA-gamma- 
Retheta) 


Page 22 


LES Models 


¢ Specifications 
¢ Large-eddy simulation (LES) based on theory of self 
similarity 
* Solves small eddies with a sub-grid-scale model and 


large ones based on geometry calculation 
¢ Separates velocity field into resolved (large eddies) 


and sub-grid parts (small eddies) 


¢ Applications 


¢ Thermal fatigue, vibration, buoyant flows (ship design) 


¢ Limitations 
¢ Difficulties in near-wall regions 
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LES Models 
Altair Autodesk 
mnouels AcuSolve CFD 
Smagorinsky X xX 


sub-grid scale 


Dynamic 
Smagorinsky xX 
sub-grid scale 
Differential SGS 
Stress Equation 
One Equation 
Eddy-Viscosity 
Dynamic One 


Equation Eddy- 
Viscosity 
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DES, DDES, IDDES Models 


¢ Specifications 


D It 
ESI Group / oe 
Systémes 


Pa patente SIMULIA 
pen ABAQUS/CED 


X X 
X xX 

xX 
X X 


* Detached eddy simulation (DES) RANS near walls and 


LES elsewhere, address LES near-wall issues 
* Delayed detached eddy simulation (DDES) 
¢ Avoid modeled-stress depletion and grid-induced 


separation 


¢ Improve delayed detached eddy simulation (IDDES) 
¢ Similar to DDES and solves log-layer mismatch 


¢ Applications 


¢ External aerodynamics, aeroacoustics, bulk or wall 


turbulence 


¢ Limitations 


¢ Lots of programming 
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DES, DDES, IDDES Models 


Models Altair Autodesk COMSOL 
AcuSolve CFD Multiphysics 
DES Spalart- 
Allmaras xX 
Detached Eddy 
DDES SST k-o 
Delayed xX xX 
Detached Eddy 
DDES Spalart- 
Allmaras Delayed XxX 
Detached Eddy 
Shielded 
Detached Eddy 
Simulation (SDES) 
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DES, DDES, IDDES Models 


Models Altair Autodesk COMSOL 
AcuSolve CFD Multiphysics 
Stressed Blended 
Eddy Simulation 
(SBES) 
IDDES SST K-o 
Improved 
Delayed xX 
Detached Eddy 
Simulation 
IDDES Spalart- 


ESI Group / 


ANSYS_ CED Direct 


XxX 


OpenFOAM 


XxX 


ESI Group / 


ANSYS_ CFD Direct 


xX 


OpenFOAM 


Dassault 

Systémes 

SIMULIA 
ABAQUS/CFD 


Dassault 

Systémes 

SIMULIA 
ABAQUS/CFD 


TAN POOERE xX xX 
Delayed 
Detached Eddy 
Simulation 
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Other Models: Direct Numerical Simulation (DNS) 


¢ Specifications 
¢ Navier Stokes equations solved numerically without 


turbulence model 
¢ Turbulence spatial and temporal scales resolved in 


mesh 


¢ Applications 
¢ Fundamental research and development of other 


turbulence models for practical applications 


¢ Limitations 
¢ Exceptionally high computational cost, low Reynolds 


numbers only, large-scale resolution 
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Other Models: Scale Adaptive Simulation (SAS) 


¢ Specifications 
¢ Similar results as LES in unsteady flow regions 
¢ Similar results as RANS models in stable flow regions 


¢ Applications 


* Separation zones, aerospace, channel pipes 


¢ Limitations 
¢ Underlying turbulence model must be unsteady, cannot be 
made unsteady with grid refinement 
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Other Models 
Dassault 
ESI G / 
Model Altair Autodesk COMSOL ANSYS CFD ne i Systémes 
aoa AcuSolve CFD Multiphysics e a SIMULIA 
pen ABAQUS/CFD 
Direct Numerical 
Simulation (DNS) x x 
Scale Adaptive xX 
Simulation (SAS) 
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Glossary of CFD Terms 


Page 31 


Glossary of Terms 


CAE: computer-aided engineering 

CFD: computational fluid dynamics 
DES: detached eddy simulation 

DDES: delayed detached eddy simulation 


DNS: direct numerical simulation 


€: dissipation rate of turbulent kinetic energy 


IDDES: improved delayed detached eddy simulation 


k: the kinetic energy per unit mass of turbulent fluctuations 
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Glossary of Terms 


LES: large-eddy simulation 


L-Vel: length-velocities 


w: specific rate of dissipation of kinetic energy 


RANS: Reynolds-averaged Navier-Stokes equations 
SAS: scale adaptive simulation 


Y+: dimensionless wall distance 
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Learn more about simulation at 
ENGINEERING.com 


Read some simulation news at: 


http://www.engineering.com/DesignSoftware/DesignSoftwareArticles/articleType/Tag View/Tag/Simulation.aspx 


Watch some simulation webinars at: 


http://www.engineering.com/Webinars.aspx 


Read more simulation resources and white papers at: 


http://www.engineering.com/Resources.aspx?sec=6 


